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Ionized iron Ka lines in AGN X-ray spectra 
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Abstract. The Equivalent Widths (EW) of the He- and H-like iron lines produced in photoionized, circumnuclear 
matter of AGN are calculated with respect to both the reflected and the total continua. We found that the EWs 
with respect to the total continuum may be as large as a few tens of eV, making them observable in bright Seyfert Is 
by instruments on board Chandra and XMM-/V ewton. We apply our calculations to the XMM-JVewton spectrum 
of NGC 5506 and found a good agreement with the data at the expense of a modest iron overabundance. 
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1. Introduction 

The basic idea behind unification models for Seyfert galax- 
ies is that all objects are intrinsically the same and all 
observed different features are simply the result of the 
orientation of the system with respect to the line of sight. 

In this scenario, emission from circumnuclear regions 
is much easier to observe in Seyfert 2s (and in particu- 
lar in Compton-thick ones), where the nuclear radiation 
is obscured, than in unobscured Seyfert Is in which these 
components are heavily diluted by the photons from the 
nucleus. Indeed, such regions are commonly observed in 
Seyfert 2s (e.g. Matt 2001 for a review), and found to be 
of at least two types: a cold gas with high column density 
which partially or completely obscures the nuclear radia- 
tion (the 'torus') and an ionized, optically thin medium 
that scatters it into the line of sight. Both these materials 
produce a reflected continuum and narrow iron Kq lines 
(e.g. Matt et al. 1996, hereinafter MBF96). 

If these emitting regions are also present in the envi- 
ronment of Seyfert Is, as expected in unification schemes, 
we should in principle be able to observe these lines as well, 
although with a much smaller EW due to the dilution by 
the direct nuclear continuum. The limited sensitivity and 
energy resolution of past X-ray missions prevented, un- 
til a couple of years ago, the unambiguous detection of 
these features in unobscured sources. The situation has 
been radically changed with the launch of Chandra and 
XMM-7V ewton. Indeed, a narrow neutral iron line have 
been found to be fairly common in Seyfert Is spectra 
(e.g. Kaspi et al. 2001; Yaqoob et al. 2001; Weaver 2002; 
Yaqoob et al 2002 for Chandra results; Reeves et al. 2001; 
Pounds et al. 2001; Gondoin et al. 2001; Matt et al. 2001 
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for XMM-Newton findings). On the other hand, evidence 
for ionized iron Kct lines at 6.7 and 6.96 keV from distant 
matter have been reported so far in only one unobscured 
(at the lines energies) Seyfert, namely NGC 5506 (Matt 
et al. 2001). 

In this paper, we calculate the expected equivalent 
widths for the two ionized lines with respect to both re- 
flected and total continua. Details of the calculations and 
results will be presented in Sects. 2 and 3, respectively, 
and the case of NGC 5506 discussed in Sect. 4. 

2. The model 

The line emitting material is assumed to be optically thin 
to Thomson scattering and in photoionization equilib- 
rium, as observed e.g. in the 'hot reflector' in NGC 1068 
(see Guainazzi et al. 1999; Bianchi et al. 2001). In these 
conditions, both recombination and resonantly scattered 
Ka lines from H- and He-like iron ions can be produced. 

Resonantly scattered Kcv lines consists of radiative 
de-excitation following excitation from continuum pho- 
tons. The radiative channel has a probability of occur- 
rence which, for Fe xxv and Fe XXVI, is unity as Auger 
de-excitation is of course impossible. The angular depen- 
dence of the emission depends on the transition involved 
(Chandrasekhar 1960). In the column density intervals ex- 
plored in this paper, resonant lines are optically thick. 
Therefore, the first resonant scattering occurs preferen- 
tially close to the illuminated surface (this effect of course 
becoming more and more important as the column den- 
sity of the matter increases). As there are not competing 
processes apart from Compton scattering (which has a 
much smaller cross section), the line photons eventually 
escape from the matter, but with a larger probability to 
escape from the illuminated surface (reflection case, see 
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Fig. 1) rather than from the opposite boundary (trans- 
mission case). Line transfer in the matter is treated here 
by means of Monte Carlo simulations. 

Recombination Ka lines occur after a photoionization 
of a K shell electron. Similarly to the case of fluorescence, 
it is customary to describe the process by introducing an 
effective fluorescent yield Y, which is the probability that 
the recombination cascade does include the Ka transition 
(Krolik & Kallman 1987; MBF96). Recombination occurs 
almost homogeneously in the matter, as it follows pho- 
toabsorption, which is optically thin in the column density 
ranges considered here. Of course, Ka line photons pro- 
duced by recombination may in turn be resonantly scat- 
tered before escaping, but because of the uniform source 
function there is no preferred escaping surface. This pro- 
cess can be therefore treated in the optically thin approx- 
imation, thus analytically. 

Both processes and their treatment are described in 
detail in MBF96, to which the reader is deferred. 

Ionization fractions for Fe xxv and Fe xxvi have been 
calculated with the well known public photoionization 
code CLOUDyQ (v94.00; Ferland 2000), as a function of 
the ionization parameter U , defined as: 



U 



X 



47rr 2 cn e 



(i) 



where c is the speed of light, r the distance of the gas 
from the illuminating source, n e its density and vr the 
frequency corresponding to 1 Rydberg. This, however, is 
not the best choice when investigating highly ionized iron 
atoms, as the result depends dramatically on the spectral 
shape of the photoionizing continuum. Assuming a power 
law shape, a small change in the spectral index (for the 
same U) may change significantly the number of X-ray 
photons (which are the most relevant in this case). We 
decided to use instead a "X-ray" ionization parameter, 
defined over the 2-10 keV energy range: 
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where the integration limits are the frequencies correspon- 
dent to 2 and 10 keV. The conversion factor between 
the two ionization parameters is clearly a function of the 
power index of the continuum: 



JVr hv 



(3) 



where 7 is the photon index of the powerlaw and Er = 
13.6 x 10~ 3 keV. Adopting this new definition of the ion- 
ization parameter, the results are almost independent of 
7: the residual effects will be presented at the end of the 
next section, together with those due to changing the tem- 
perature and density of the gas. 
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Fig. 1. The adopted geometry. 

The sketch in Fig|l] shows the geometry adopted in 
our calculations. We assumed Ar=r and isotropic illumi- 
nation within a 30° cone. As said above, the equivalent 
widths of the resonant scattered lines are different when 
transmitted or reflected spectra are considered. Therefore, 
we will present plots for both cases. 

3. Results 

In Figs. ||, ||, d and || the equivalent widths of Fe xxv 
and Fe xxvi Ka lines (recombination plus resonantly scat- 
tered; hereinafter we will always consider the sum of the 
two contributions) against the transmitted (or reflected) 
continuum only are plotted as a function of log U x and for 
different values of the column densities of the line emitting 
gas. The results are obtained with 7 = 2.0, T — 10 6 K, 
n e = 10 6 cmT 2 . Plots in Figs. || and || were obtained com- 
bining the results in figure 4, 5 and 6 of MBF96 (where 
EWs of recombination and resonant lines where presented 
assuming ionization fractions equal to 1 for the interested 
ion) with the ionization structures calculated by CLOUDY 
as a function of U x . The iron abundance is assumed equal 
to the solar one (as for Morrison & McCammon 1983). As 
described in MBF96, for small column densities the main 
contribution comes from the resonant lines, while at large 
column densities recombination dominates, because pho- 
toionization is still optically thin while resonant scatter- 
ing is thick, and the line production for the latter process 
consequently saturates. Plots in Figs. || and || have been 
obtained re-running the Monte Carlo code described in 
MBF96 for a reflecting, rather than transmitting, mate- 
rial. Because the two cases differs only for the contribution 
of the resonant lines, the main differences occur for small 
column densities. 

In Figs. H and the EWs, obtained for the reflecting 
matter, are instead calculated against the total contin- 
uum. (The corresponding plots for the transmitted flux 
are not shown because differences are appreciable only for 
small column densities, where the EWs are very small.) 
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Fig. 2. Fe xxv EWs against transmitted continuum only 
as a function of log U x and column density. 
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Fig. 4. Fe xxv EWs against reflected continuum only as 
a function of log U x and column density. 

The relation between the EWs with respect to reflected 
and total continua is easy to calculate. In the Compton- 
thin regime, the reflected flux is given by: 

F r ~f-T-F t (4) 

where F t is the illuminating flux and / a geometrical 
factor (basically the solid angle subtended by the illumi- 
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Fig. 3. Fe xxvi EWs against transmitted continuum only 
as a function of log U x and column density. 
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Fig. 5. Fe xxvi EWs against reflected continuum only as 
a function of log U x and column density. 

nated matter); t = 1 ^iq-m is the Thomson optical depth. 
Therefore, the equivalent widths are: 

EW t = EW r ■ ~EW r ~2iEW r -f-T (5) 

Figs. ^ and have been obtained assuming /=1. 

Finally, we illustrate the effects of changing the value 
of 7, T and n e . Fig. [s] shows the equivalent widths of 
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Fig. 6. Reflected Fe xxv EWs against total continuum 
as a function of \ogU x and column density (f=l: see text 
for details). 

Fe xxv lines against the total continuum for a column 
density of 10 23 cm -2 and different values of the photon 
index of the incident continuum. Clearly, the choice of U x 
instead of U keeps the differences small: the maximum of 
the function occurs almost at the same value of U x for 
the various indices, while it would have been significantly 
shifted if plotted against U. There is however a difference 
in the value of the maximum, which increases with 7. 

The opposite occurs for Fe xxvi (Fig. ||), with a shift 
in the value of U x for which the function has the maximum 
value, but not appreciable differences in the values of this 
maximum. 

The effect of changing the gas temperature is shown in 
Figs. [To] and [fl[ The overall result is, as expected, that the 
colder is the gas the higher is the ionization parameter re- 
quired to reach the peak in the ionization fraction of both 
the iron ions. On the other hand, there are no significant 
changes in the maximum values of the equivalent widths 
of the lines. 

A change of the gas density over a wide range (10 4 — 10 8 
cm~ 3 ) does not instead produce any appreciable variations 
in our results. This is expected, because the ionization pa- 
rameter is a quantity that describes the ionization struc- 
ture of a gas under a combination of density, distance and 
luminosity of the incident radiation. Changing one of these 
quantities keeping the ionization parameter fixed ends up 
only in a rescale of the other ones, without changing the 
ionization structure of the gas surface. 



Fig. 7. Reflected Fe xxvi EWs against total continuum 
as a function of log U x and column density (f=l: see text 
for details). 

4. An example: NGC 5506 

Let us now compare our calculations with the only case in 
which the presence of ionized iron Ka lines from circum- 
nuclear matter have been claimed in an unobscured (at 
the iron line energies) source, namely NGC 5506 (Matt et 
al. 2001). 

NGC 5506 is a bright Narrow Emission Line Galaxy, 
optically classified as an intermediate Seyfert with a highly 
reddened BLR (Vcron-Cctty & Veron 2001). Several X- 
ray observations clearly indicate the presence of a hidden 
Seyfert 1 nucleus (Matt et al. 2001 and references therein). 
Its iron Ka line has been claimed to be complex by ASCA 
(Wang et al. 1999) and RXTE (Lamer et al. 2000) ob- 
servations. BeppoSAX found an iron line centroid energy 
definitely higher than 6.4 keV (Perola et al. 2002). 

XMM-Newton data definitely confirm that the iron 
line is complex, with at least two components clearly 
present: one narrow, corresponding to neutral iron, the 
other broad and consistent with ionized iron. 

Even if a fit of the ionized, broad component with a 
disc relativistic line is as good as that with a blend of 
He- and H-like unresolved lines, the second solution is 
preferable from a physical point of view, as argued by 
Matt et al. (2001). If this is, indeed, the case, these two 
lines are likely to be produced in a photoionized matter as 
the one described in this paper. The observed equivalent 
widths are 40 ± 16 eV and 32 ± 15 eV for Fe xxv and 
Fe xxvi, respectively. 

Moreover, NGC 5506 presents a soft excess whose nor- 
malization is about 2% of the primary component: recall- 
ing equation (^) this could be the result of reflection from 
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Fig. 8. Fe xxv EWs against total continuum for different 
values of 7. Column density is 10 23 cm~ 2 . 
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Fig. 10. Fe xxv EWs against total continuum for differ- 
ent values of T. Column density is 10 23 cm~ 2 . 

a gas with column densities of 3 x 10 22 // cm~ 2 . Let us 
now assume that the two ionized lines and the soft X-ray 
emission come from one and the same matter. Let us also 
assume f=l (a smaller value of / would implies a larger 
value of the column density to produce the observed flux 
of the reflected continuum: the two parameters have the 
same effect on the line EWs, and in the first approxima- 
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Fig. 9. Fe xxvi EWs against total continuum for different 
values of 7. Column density is 10 23 cm~ 2 . 
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Fig. 11. Fe xxvi EWs against total continuum for differ- 
ent values of T. Column density is 10 23 cm~ 2 . 

tion we can assume that the EWs are almost independent 
of their ratio). The observed EWs can be reproduced pro- 
vided that \og{U x ) ~ 0.5 (T = 10 6 K) and the iron is 
overabundant by a factor 2-^3. (A different choice of the 
gas temperature would require a different value of log(Z7 x ), 
but not of the iron abundance, see Figs. 10 and 11). This 
solution is of course valid in the stationary case. If the 
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source flux was larger in the past than during the XMM- 
Newton observation, then the requirement on the iron 
overabundant will be correspondingly reduced; it is worth 
noting that the flux during the XMM-Newton observa- 
tion was about 30% lower than the average flux during 
the longer contemporaneous BeppoSAX observation. 

Finally, combining the derived value of the ionization 
parameter with the luminosity of the object (£2-10 = 
8.5 x 10 42 erg cm 2 s -1 ), its spectral shape (T ~ 2) and 
the inferred column density of the gas, we can write two 
equations linking the above, observed parameters with the 
density n e , the distance r (from the nucleus) and the depth 
Ar of the gas (see Fig. : 



10 
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hv 
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where r = Ar in our model (see Sect. |J). If we use the 
value of U x found before to reproduce the observed lines 
for T = 10 6 K, we get the following estimates for r and 



2.5 x 10 5 // 



-3 



0.04-/ pc 



(7) 



These values change if we use the ionization parame- 
ters which best reproduce the observed EWs in the case 
of T = 10 5 K and T — 10 7 K (see Figs. and [□]), being 
respectively: [n e ~ 8 x 10 5 // cm~ 3 ; r ~ 0.01 • / pc] and 
[n e ~ 8 x 10 4 // cm" 3 ; r ~ 0.1 • / pc]. 

This distance scale is slightly larger than that typical 
of the BLR (e.g. Peterson 1997). On the other hand, it is 
definitely smaller than the inner radius of the absorbing 
torus, which is of the order of parsecs (e.g. Bianchi ct al. 
2001). Interestingly, our estimated distance is consistent 
with the (poorly constrained) distance scale of warm ab- 
sorbers (e.g. Otani et al. 1996; Netzer et al. 2002), but it is 
worth remarking that the ionization parameter required in 
NGC 5506 is clearly higher than that of warm absorbers. 

5. Conclusions 

We have calculated the equivalent widths of H- and He- 
like iron Ka lines from photoionized, Compton— thin mat- 
ter in Seyfert galaxies, using results from MBF96 and ion- 
ization fractions estimated with CLOUDY. We have applied 
our calculations to NGC 5506, whose XMM-Newton spec- 
trum shows a complex iron Ka line system, possibly com- 
posed by a neutral line and a blend of two ionized lines 
as those produced in our model. We found that, if this 
interpretation is correct, the ionized lines observed in this 
source can be reproduced assuming log U x ~ 0.5 (if 
T = 10 6 K) and an iron overabundance of 2 3. 

At a given gas temperature, the range of ionization pa- 
rameters which produces EWs as observed in NGC 5506 is 
rather narrow; moreover, NGC 5506 is one of the bright- 
est AGN in the X-ray sky. We then expect that, unless 
for some unknown reason the ionization parameter tends 
to assume values in this narrow range, in most cases these 



lines will be below detectability even with XMM-Newton 
and Chandra, and therefore that NGC 5506 may be the 
exception rather than the rule. It is interesting to note 
that very similar iron line complexes have been observed 
so far by XMM-Newton in only two other sources, Mrk 205 
(Reeves et al. 2001) and Mrk 509 (Pounds et al. 2001) 
where, as discussed in those papers, the origin of the ion- 
ized lines is likely to be the accretion disc. 
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